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MEASUREMENTS  OF  LASER  EXTINCTION  IN  ICE  FOG  FOR 
DESIGN  OF  SEV  PILOTAGE  SYSTEM 


by 

Richard  Munis  and  Allan  Delaney 

INTRODUCTION 

The  operation  of  any  vehicle  during  conditions  of  low  visibility  ic-'uires  that  the  operator  be 
able  to  “see”  through  various  weather  conditions,  including  fog,  rain  and  snow.  The  operation  of 
the  vehicle  on  arctic  terrain  in  addition  requires  that  the  operator  be  prepared  to  avoid  ridges  on 
sea  ice.  This  combination  of  problems  must  be  considered  in  the  design  of  an  obstacle  avoidance 
system  for  a  surface  effect  vehicle  (SEV). 

The  Advanced  Research  Projects  Agency  (ARPA)  has  initiated  a  program  on  electromagnetic 
propagation  to  study  the  effects  on  vehicle  operation  of  snow  and  ice  covered  terrain,  and  ice  and 
water  fog  at  visible,  infrared,  millimeter,  and  microwave  wavelengths.  As  part  of  this  program,  a 
report  has  been  prepared  on  the  radar  cross-section  measurements  of  snow  and  ice,*  The  purf  so  of 
the  present  report  is  to  present  data  on  laser  extinction  measurements  in  ice  fog  at  wavelengths  of 
0.6328,  1.15  and  3.39  g. 

LITERATURE  REVIEW 

x 

A  great  deal  of  effort  has  been  spent  in  studying  the  mechanism  of  ice  fog  formation.  In  1954, 
Thuman  and  Robinson5  made  the  first  observations  of  the  size  distribution  of  ice  fog  crystals  near 
Fairbanks,  Alaska.  In  1964,  Kumai5  prepared  a  report  on  the  concentration  and  size  distribution 
measurements  of  ice  fog  also  in  the  vicinity  of  Fairbanks.  Most  recently,  Benson*  and  Weller7  sum¬ 
marized  the  results  of  these  investigations  of  the  Fairbanks  ice  fog  and  Ohtake*  attempted  to  con¬ 
firm  the  mechanism  of  ice  fog  formation  that  he  and  other  investigators  had  previously  postulated. 

Arv'.  ding  to  Ohtake,*  the  primary  sources  of  ice  fog  in  the  Fairbanks  area  were  open  water,  and 
automobile  and  heating  plant  exhaust.  One  might  at  first  question  whether  such  sources  of  ice  fog 
would  be  significant  to  a  SEV  operating  in  arctic  regions,  well  away  from  the  influence  of  city 
environments;  however,  it  is  quite  possible  that  ice  fog  would  form  on  nuclei  from  the  SEV  exhaust. 
Thus,  any  suitably  designed  obstacle  avoidance  system  for  the  SEV  should  be  provided  with  the  r 
capability  of  detecting  obstacles  through  low-visibility  conditions  generated  by  the  vehicle  itself. 

Although  much  theoretical  and  experimental  research  has  been  done  on  the  scattering  and 
absorption  properties  of  warm  fog,  very  little  work  has  been  done  on  the  extinction  properties  of  ice 
fog.  One  significant  reason  for  the  lack  of  data  may  be  the  difficulty  of  producing  ice  fog  in  the 
laboratory.  Another  may  be  the  problem  of  obtaining  f  ield  data  at  temperatures  as  low  as  -40°C. 
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Regardless  of  the  reason  the  present  investigation  is  not  based  on  the  extensive  experience  of 
other  workers.  The  various  problems  confronting  the  researcher  who  attempts  to  make  meaningful 
extinction  measurements  in  ice  fog  and  to  relate  them  to  existing  theory  are  discussed  in  Appendix 
A. 


EXPERIMENTAL  PROCEDURE 

The  extinction  measurements  at  0.6328,  1.15  and  3.39  g  were  made  in  a  4-m-long  fog  chamber 
located  in  a  cold  room  whose  temperature  could  be  lowered  to  approximately  -36°C.  Temperatures 
below  — I0°C  were  achieved  inside  the  fog  chamber  by  cooling  coils  located  inside  the  chamber. 
Figure  1  shows  an  exterior  view  of  this  chamber  and  Figure  2  shows  the  optical  path  down  which 
the  laser  beam  was  propagated  toward  the  detector.  Figure  3  shows  the  interior  of  the  room  in 
which  the  helium-neon  laser  and  the  electronic  processing  equipment  were  located.  The  tempera¬ 
ture  of  this  room  was  kept  at  approximately  20°C. 

The  extinction  measurements  were  made  in  the  following  manner.  The  modulated  laser  beam 
was  sampled  at  the  opposite  end  of  the  4-m  fog  chamber  by  a  photodiode  and  the  signal  obtained  was 
processed  by  a  lock-in  amplifier  before  being  recorded.  The  strip  chart  recorder  was  allowed  to  run 
long  enough  to  establish  a  stable  signal  level.  This  signal  level  was  then  used  as  a  measure  of  the 
incident  radiation  impinging  upon  the  fog. 

Following  this,  a  fog  was  produced  in  the  chamber  by  using  a  technique  that  is  discussed  in 
Appendix  B.  As  the  attenuation  of  the  laser  beam  was  continuously  recorded,  fog  particles  were 
sampled  with  an  impactor  (Fig.  4).  The  fog  concentration  varied  over  a  rather  wide  range  while  the 
beam  attenuation  was  being  measured.  Immediately  after  the  fog  particles  were  sampled,  the  glass 
slide  was  removed  from  the  impactor,  placed  under  a  microscope  and  photographed.  Figures  Cla- 
C7a,  Appendix  C,  show  ice  fog  samples  taken  during  these  experiments.  Photographs  of  these  sam¬ 
ples  were  later  enlarged  to  obtain  particle  spectra.  Figures  Clb-C7b  show  the  corresponding  plots 
of  diameter  d  vs  number  of  particles  N  for  each  of  the  ice  fog  samples. 


Figure  I.  Exterior  view  of  fog  chamber,  fee  fog  sampler  a/e  obtained  hv 
inserting  the  impactor  into  one  of  the  circular  ports. 


propagation.  Over- 
to  the  chamber. 


laser  and  electronic  proc- 


ft 


Figure  4.  Impactor  with  which  log  samples  were  obtained.  As  a  measured 
volume  ol  air  is  drawn  through  on  orifice  by  pulling  back  the  syringe,  ice 
fog  particles  contained  in  the  air  sample  are  deposited  in  a  film  of  silicone 
oil  on  a  microscope  slide  located  underneath  the  orifice. 


EXPERIMENTAL  AND  THEORETICAL  EQUATIONS 


The  extinction  coefficient  at  each  wavelength  a  for  each  particle  concentration  N  may  be  cal¬ 
culated  according  to 


T(A,  N)  =  e~a(A'  V>f 

where 

T(A,  N)  =  measured  transmission  (%) 
l  =  optical  path  length  (m) 
a(A,  N)  =  extinction  coefficient  (m-1). 

The  extinction  coefficient  for  a  single  particle  may  be  written  as.6 

u  =  S  +  A 


(1) 


(2) 


where 

S  =  scattering  coefficient  (m-1) 

A  =  absorption  coefficient  (m-1). 

At  a  wavelength  of  0.6328  p  absorption  is  insignificant  because  most  of  the  radiation  not  impinging 
upon  the  detector  is  lost  only  to  angular  scattering  out  of  the  primary  beam.  However,  at  3.39  \i  both 
terms  in  eq  2  contribute  to  the  loss  of  radiation  as  the  beam  traverses  the  fog  chamber.  Thus,  at 
this  wavelength  the  complex  part  of  the  index  of  refraction  becomes  a  significant  factor  in  calculating 
the  theoretical  extinction  coefficient. 

The  magnitude  of  the  theoretical  extinction  coefficient  may  be  calculated  from  the  particle  size 
distribution  and  the  Mie  theory  using  the  relation6 

a 

max  p 

y=  X  N(a)  Aaira^Qext(m,x) 
amin 


(3) 
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where 

N( a)  =  particle  concentration  per  Aa  radius  interval 

ra2  =  geometrical  cross-section  area  of  a  single  particle 
<?ext  =  van  de  Hulst’s  dimensionless  efficiency  factor  fa  total  extinction* 

m  =  complex  index  of  refraction 
x  =  particle  size  parameter  2na/k. 

Equation  3  is  rigorously  valid  only  fa  single  scattering  caused  by  spherical  particles. 


RESULTS  AND  DISCUSSION 

Table  I  compares  the  theoretical  and  experimental  extinction  coefficients  erf  ice  fog  at 
0.6328  /x  computed  from  eq  1  and  3.  Also  included  are  the  total  particle  concentration  of  each  sample 
and  the  measured  value  of  transmission.  The  last  column  of  this  table  gives  the  aitenuation  in  terms 
of  db/m,  a  value  that  is  probably  mote  useful  to  the  systems  designer. 


Table  I.  Experimental  and  theoretical  extinction  coefficients 
of  ice  fog  at  0.6328-/1  and  -36° C. 


Panicle 

concentration 

(N/cms) 

Experimental 

extinction 

coefficient 

(m~l) 

Mie 

extinction 

coefficient 

On’1) 

Transmission 

(%) 

Attenuation 

(db/m) 

66 

0.0320 

0.0212 

88 

-0.0138 

20 

0.0496 

0.0095 

82 

-0.0215 

15 

0.0621 

0.0126 

78 

-0.0268 

108 

0. 1236 

0.0268 

61 

-0.0534 

Figure  5  shows  the  behavior  of  the  experimental  extinction  coefficient  at  0.6328  y  as  the 
particle  concentration  increases  and  the  predicted  behavior  of  the  extinction  coefficient  according 
to  Mie  theory.  Although  there  is  some  disagreement  between  theay  and  experiment,  both  indicate 
the  same  behavia  as  particle  concentration  increases;  that  is,  after  a  certain  particle  concentra¬ 
tion  has  been  attained  asymptotically  no  further  increase  of  the  extinction  takes  place.  During  this 
set  ot  measurements  the  authors  experienced  significant  difficulty  in  generating  a  repeatable  ice  fog; 
therefore,  these  data  had  to  be  assembled  from  two  different  fog  samples.  Fig.  Cla  and  Clb  show 
the  fog  sample  and  the  corresponding  particle  spectra  for  one  of  these  samples. 

Tables  II  and  III  give  the  experimental  and  theaetical  extinction  coefficients  of  ice  fog  at 
1.15  y  and  3.39  y  respectively  for  a  chamber  temperature  ot  -43°C.  Figures  6  and  7  show  the  be¬ 
havia  of  the  two  coefficients  as  the  particle  concentration  increases.  At  1.15  y  the  agreement 
between  the  two  is  relatively  good.  However,  the  most  striking  aspect  of  this  comparison  is  that 
both  experiment  and  theay  indicate  a  change  of  the  rate  of  extinction  at  about  90  particles/cm2. 

At  this  value  a  new  rate  is  established  and  both  experiment  and  theory  show  approximately  the  same 
general  trend  for  the  available  data.  By  this  time  the  prohlem  of  generating  a  repeatable  fog  had  been 
overcome;  thus  these  data  are  from  samples  taken  from  one  fog  as  it  dissipated. 


Figure  5.  Experimental  and  theoretical  extinc-  Figure  6.  Experimental  and  theoretical  extinc 

don  coefficients  as  particle  concentration  don  coefficients  as  particle  concentration  in¬ 
increases  at  0.6328  p.  creases  at  1.15  p. 


Figure  7.  Experimental  and  theoretical 
extinction  coefficients  as  particle  con¬ 
centration  increases  at  3.39  p. 
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Table  n.  Experimental  and  theoretical  extinction  coefficients 
of  ice  fog  at  1. 15  p  and  -43°C. 


Particle 
concentration 
( N/cm ') 

Experimental 

extinction 

coeflicient 

(m'1) 

Hie 

extinction 

coefficient 

On'1) 

Transmission 

(%) 

Attenuation 

(db/m) 

70 

0.0128 

0  0214 

95 

-0.0055 

91 

0.0234 

0.0401 

91 

-0.0102 

159 

0.032 

0.0649 

88 

-0.0138 

Table  HI.  Experimental  and  theoretical  extinction  coefficients 
of  ice  fog  at  3.39  p  and  -43°  C. 


Particle 

concentration 

(N/cm') 

Experimental 

extinction 

coefficient 

(m'1) 

Mie 

extinction 

coefficient 

(m‘l) 

Transmission 

(%) 

Attenuation 

(db/m) 

54 

0.0075 

0.0182 

97 

-0.00324 

165 

0.0128 

0.0641 

95 

-0.0055 

101 

0.0128 

0.0432 

95 

-0.0055 

238 

0.0207 

0. 1024 

92 

-0.00895 

155 

0.0207 

0.0628 

92 

-0.00895 

209 

0.0348 

0.0731 

87 

-0.0150 

469 

0.0557 

0.2027 

80 

-0.024 

Figure  8  shows  a  composite  plot  of  the  particle  spectra  at  sampling  times  t0,  tj  and  t2.  At 
t0  the  peak  of  the  distribution  was  located  at  15.2  p,  at  tj  it  had  shifted  slightly  to  14.4  p,  and  at 
t2  a  significant  shift  located  it  at  11.6  p.  Figure  8  also  shows  that  at  t0  the  particle  concentration 
was  165/cm3,  but  when  it  had  reached  91/cm3  the  spectrum  shifted  such  that  there  were  fewer 
particles  and  the  mean  particle  size  was  getting  smaller.  At  this  point  the  rate  of  extinction  suddenly 
increased  and  the  mean  diameter  of  the  particle  shifted  again  to  an  even  smaller  size.  It  is  not 
possible  at  this  time  to  determine  whether  the  extinction  coefficient  at  70/cm3  was  smaller  than  that 
at  165/cm3  because  the  number  of  particles  was  smaller  oi  because  the  mean  diameters  of  the 
particles  were  smaller. 

Figure  9  shows  a  composite  plot  of  the  particle  spectra  taken  during  the  experimental  work  at 
3.39  p.  Here  a  different  distribution-time  relationship  evolved;  that  is,  the  particle  spectrum  first 
shifted  to  a  larger  mean  size  before  it  shifted  back  again  to  the  original  distribution,  where  it  remained 
unchanged  until  the  fog  dissipated.  Figure  10  shows  a  plot  of  concentration  versus  extinction  coeffi¬ 
cient  for  the  data  samples  obtained  during  the  same  fog.  Unfortunately  agreement  between  theory 
and  experiment  was  not  as  remarkable  as  at  1.15  p. 

Figure  11  plots  the  experimental  extinction  coefficient  in  units  of  db/m  versus  particle  con¬ 
centration  at  the  three  laser  wavelengths.  Again,  it  must  be  emphasized  that  more  than  one  particle 
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size  distribution  is  represented  by  these  data.  However,  the  important  feature  here  is  the  comparison 
at  these  three  wavelengths  of  the  attenuation  properties  of  ice  fog  of  varying  concentration.  From 
Figure  11  we  may  deduce  that  at  0.6328  \l  only  a  few  particles  per  cm  cause  an  appreciable  loss  of 
energy  in  the  form  of  scattering.  However,  at  the  two  infrared  wavelengths  (1.15  p  and  3.39  p)  the 
same  concentration  causes  only  a  negligible  loss  of  energy.  Further  inspection  shows  that  in 
extremely  dense  fogs  loss  of  energy  becomes  appreciable  at  3.39  whereas  at  1.15  /*  it  appears 
from  the  available  data  that  the  extinction  coefficient  doss  not  rise  as  rapidly  as  it  does  at  3.39  fi. 
This  seems  to  agree  with  the  fact  that  the  imaginary  part  of  the  index  of  refraction  of  ice  is  about 
104  higher  at  3.39  p  than  at  1.15  p. 


CONCLUSIONS 

Experimental  and  theoretical  data  have  been  presented  for  three  laser  wavelengths:  0.6328,  1.15  and 
3.39  fi.  The  authors  experienced  a  good  deal  of  difficulty  at  the  beginning  of  the  program  in  gen¬ 
erating  and  obtaining  good  ice  fog  samples.  However,  these  problems  have  been  largely  overcome 
and  it  is  now  possible  to  generate  repeatable  fogs.  The  results  of  this  study  have  shown  that  there 
is  generally  fairly  good  agreement  between  the  experimentally  derived  extinction  coefficients  and 
those  predicted  by  the  Mie  theory.  More  remarkable  is  the  fact  that  at  1.15  p  experiment  and  theory 
show  excellent  agreement  between  the  behavior  of  the  extinction  coefficient  and  particle  concentra¬ 
tion.  However,  more  data  are  needed  at  both  lower  and  higher  particle  concentrations  to  substantiate 
the  agreement  between  theory  and  experiment. 
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APPENDIX  A:  SOME  PROBLEMS  IN  THE  EXPERIMENTAL  DETERMINATION 
OF  THE  EXTINCTION  COEFFICIENT  OF  ICE  FOG 

Recently  the  difficulty  of  reliably  characterizing  an  ice  fog  distribution  by  particle  size  and 
number  has  been  markedly  reduced.  The  viscosity  of  oil  used  to  collect  ice  particles  is  very 
important  in  determining  any  distribution  because  if  the  oil  is  too  thin  and  large  particles  are 
sampled  it  cannot  support  the  momentum  of  the  particles  as  they  impact  on  the  slide  and  they  there¬ 
fore  crash  into  the  slide  and  disintegrate. 

However,  even  if  proper  oil  viscosity  is  selected  success  cannot  be  guaranteed  because  the 
problem  01  depositing  the  oil  on  the  slide  still  remains.  If  care  is  not  taken  to  restrict  the  oil  from 
spreading  too  much  after  it  is  deposited,  the  particles  still  tend  to  disintegrate  because  of  the 
"thinness”  of  the  oil  deposit  rather  than  the  viscosity.  The  authors  conducted  a  series  of  tests  in 
which  the  same  viscosity  of  oil  was  deposited  in  a  slightly  different  manner.  If  extreme  care  was 
not  taken  a  variation  of  up  to  400“'  could  be  expected  in  the  total  number  of  particles  collected  on 
the  slide.  Even  with  extreme  care  variations  of  between  10-50°.  were  observed.  Clearly  this 
problem  must  be  given  more  attention  if  better  agreement  is  to  be  obtained  between  theory  and  experi¬ 
ment. 

Another  problem  has  recently  come  to  the  authors’  attention,  if  a  rather  small  aperture  is  used 
on  the  impactor  to  admit  the  fog  particles  and  if  the  syringe  is  spring-loaded,  the  velocity  with 
which  the  air  enters  the  sample  chamber  is  so  high  that  it  causes  the  oil  to  separate,  thus  allowing 
the  particles  to  crash  into  the  glass  slide  and  disintegrate.  However,  releasing  the  syringe  slowly 
prevents  this  problem  and  makes  it  possible  to  obtain  rather  consistent  results. 
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APPENDIX  B:  ICE  FOG  PRODUCTION  TECHNIQUE 


The  fog  production  technique  uses  a  pressure  cooker  to  produce  steam  which  is  then  sprayed 
directly  into  the  fog  chamber  until  the  chamber  is  saturated.  Ice  crystals  begin  to  form  almost 
as  soon  as  the  steam  starts  to  cool  down  to  the  chamber  temperature.  At  this  point  the  fog  is 
“moved”  around  in  the  chamber  until  it  is  thoroughly  mixed.  The  reason  for  this  is  that  since  the 
steam  is  admitted  into  the  chamber  at  only  one  location  the  laser  propagation  path  tends  to  have  a 
more  random  fog  particle  distribution  along  the  complete  length  of  the  path  than  if  the  fog  were 
permitted  to  spread  from  its  initial  position.  After  the  fog  has  been  stirred  3bout  it  is  then  per¬ 
mitted  to  dissipate  without  any  further  interferences.  As  it  dissipates  samples  are  taken  and 
photographed  approximately  every  2-3  min  until  it  is  completely  dissipated;  this  is  done  in  con¬ 
junction  with  the  continuous  recording  of  the  laser  beam  attenuation. 

The  authors  have  been  greatly  encouraged  by  the  close  resemblance  of  the  fog  samples  taken 
during  this  investigation  to  those  obtained  by  Ohtake  in  Fairbanks.  Alaska.4  Figure  B1  shows  the 
range  of  the  mean  diameter  dm  of  spherical  ice  particles  versus  temperature  for  samples  taken 
during  this  investigation,  and  Figure  B2  shows  the  data  of  Ohtake4  and  Thuman  and  Robinson5  taken 
in  Alaska.  A  comparison  of  these  results  indicates  that  the  technique  for  generating  ice  fog  in  the 
environmental  chamber  used  in  the  present  investigation  is  indeed  realistic. 


Figure  Bl.  Mean  diameter  of  spher-  Figure  B2.  Mean  diameter  of  spherical  ice 

ical  ice  particles  versus  tempera-  panicles  versus  temperature.* 

lure. 
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APPENDIX  C:  ICE  FOG  SAMPLES  AND  THEIR 
CORRESPONDING  PARTICLE  SPECTRA 

(N  =  number  of  particles  and  d  =  diameter  of  ice  fog  sample.) 
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Clb. 
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